GLOBAL PARAMETRICES AND DISPERSIVE ESTIMATES FOR 
VARIABLE COEFFICIENT WAVE EQUATIONS 



JASON METCALFE AND DANIEL TATARU 



Abstract. In this article we consider variable coefBcient time dependent wave equations in 
R X M" . Using phase space methods we construct outgoing parametrices and prove Strichartz 
type estimates globally in time. This is done in the context of metrics which satisfy a 
weak asymptotic flatness condition at infinity. 



Contents 

1. Introduction E 

2. Notations 11 

3. A minor simplification 12 

4. The paradifferential calculus 14 

5. Localized energy estimates 19 

6. The half wave decomposition 24 

7. Parametrices and Strichartz estimates 26 

8. Pseudodifferential operators and phase space transforms 37 

9. A long time phase space parametrix 40 

10. A perturbation of the half wave equation 42 

11. The parametrix construction 45 
References 51 



1. Introduction 
Begin with the constant coefficient wave equation in ] 

nu = {dt - A)u = 0, n(0) = uo, Mt(0) = ui. 
On one hand the energy is preserved, 

\\Vu{t)\\L2 = \\Vu{0)\\L2 
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where V stands for the space-time gradient of the solution. On the other hand there is 
pointwise decay of waves with locahzed initial data. Precisely, if we set Uq = then 



1 — n / Nil , n — 1 



(1) \\\Dx\ 2 u{t)\\L^ < t 2 llmllii 

for all initial data Ui with a dyadic frequency localization. As a consequence of this one 
obtains the Strichartz estimates, which have the form 

(2) \\\Dx\~''Vu\\lpli < ||Vmo||l2 + ||ui||i2. 

This holds for all pairs {p,p, q) satisfying the relations 2<]9<oo,2<g<oo and 

Inn 2 n — 1 n — 1 

(3) - + - = --p, - + <— — 

p q 2 P q 2 

with the exception of the forbidden endpoint (1,2, oo) in dimension n = 3. All {p,p,q) 

satisfying these relations are called, in the sequel, Strichartz pairs. If the equality holds in 

the second part of ([3]) then the corresponding pair is called a sharp Strichartz pair. 

In the sequel, we shall not explicitly deal with the case of g = oo, allowing us to freely use 
Littlewood-Paley theory. When g = oo, we only obtain estimates with L°° replaced by the 
appropriate Besov spaces. With additional work, possibly akin to the modifications 

given to obtain the X estimate in ( l25l) . we believe that the proper estimate can be recovered. 
In what follows, we shall also concentrate our efforts on the cases that the Strichartz pairs 
are sharp. This can be done without loss of generality since the remaining estimates can 
then be recovered using Sobolev embeddings. 

A straightforward consequence of ([2]) is an estimate for solutions to the inhomogeneous 
problem 

nu = f, u{0) = 0, ut{0) = 0, 

namely 

(4) |||/^.r^Vn|UpL. < ||/||lil2. 
The simplest case of (jl]) is the well-known energy estimate 

(5) \\Vu\\l^l^ < ||/|Ui^2. 

However, there is a larger family of estimates for solutions to the inhomogeneous wave 
equation where we also vary the norms in the right hand side, 

(6) \\\D,\-fVu\\L.L'> < IllZ^.rVII^Pi^^i- 

This holds for all Strichartz pairs {p,p,q), {pi,Pi,qi)- 
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Estimates of the above type were first proved in the constant coefficient case in |3], |31j . 
Further references can be found in a more recent expository article [10]. The endpoint 
estimate (p, q) = (2, ^^3^) was only recently obtained in |13j {n > A). 

In this article we are interested in the variable coefficient case of these estimates, where 
we replace □ by a second order hyperbolic operator of the form0 

P(t, X, D) = Daa'^'^it, x)D(3 + b'^it, x)D^ + c(t, x). 

where = dk/i- Here the matrix is assumed to have signature (n, 1), and the time 
slices are assumed to be space-like, i.e. < 0. Thus we consider evolutions of the form 

(7) Pu = f, uiO) = Uo, UtiO) = u,. 

Locally in time this problem is well understood. If the coefficients are smooth then para- 
metrices are obtained using Fourier integral operators, and the Strichartz estimates were 
established in [19]. Operators with C^'^ coefficients were first considered in |25|, where a 
wave packet parametrix is constructed in all dimensions and the Strichartz estimates are 
proved in low dimension n = 2,3. An alternative parametrix construction, based on the FBI 
transform, was later obtained in [32|, [33] [35|. There the Strichartz estimates are obtained 
first for C^'^ coefficients and then for V^a G L^L°°. Below this regularity threshold for the 
coefficients the full Strichartz estimates are lost (see [2S|,[2B])! and one only retains partial 
results (see [53], [53]). 

Our goal here is to study the global in time behavior, which is a considerably more difficult 
problem. The present article is inspired by an earlier article of the second author [37| which 
deals with the same issues for the corresponding Schrodinger equation. There are many 
similarities between the two problems, but also differences. In what follows we try to discuss 
both problems in parallel. 

The dynamics for high frequencies are closely related to the Hamilton flow dynamics, 
although perhaps less so than in the case of the local in time problems. 

A first phenomena that one needs to consider is that of refocusing, which in general 
precludes the dispersive estimates ([1]) even if we restrict ourselves to coefficients which are 
sufficiently small, smooth, compactly supported perturbations of the (Minkowski) identity. 
This is because even a small perturbation of the fiat metric suffices in order to refocus a 
group of Hamilton flow rays originating at the same point and thus produce caustics. 

''^Here we employ the summation convention where repeated indices are imphcitly summed. Repeated 
Greek letters a,(3, . . . are summed from to n, where Dq — Dt, and repeated Latin indices . . . are 
summed from 1 to n. 
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At the parametrix level this is reflected in the fact that a good parametrix along a ray 
which crosses through a bounded region is very difficult to construct. This is why, following 
[37], we construct an outgoing parametrix, which only requires the analysis of the outgoing 
Hamilton flow. The price we pay is that our parametrix cannot evolve only forward in time; 
instead it must have a forward and a backward component. 

In the case of the Schrodinger equation, this is seen on arbitrarily small time scales due 
to the infinite speed of propagation; for the wave equation, on the other hand, a large time 
scale is needed. 

A second feature is related to the long time behavior of the bicharacteristics. In the fiat 
case all bicharacteristics are straight so they escape to infinity both forward and backward 
in time. However, in the variable coefficient case, it is possible to have trapped rays, which 
are confined to a bounded spatial region. These correspond to singularities which are largely 
concentrated in a bounded region and may destroy not only the dispersive estimates ([1]) but 
also the Strichartz estimates in ([21). On the positive side, the nonexistence of trapped rays 
is a more stable phenomena; in particular, it cannot happen for small perturbations. Again, 
this obstruction is seen even on short time scales for the Schrodinger equation, but only on 
large time scales for the wave equation. 

The local in time problem for the Schrodinger equation has been previously considered by 
other authors. Stafillani and Tataru [29] study compactly supported perturbations of the 
fiat metric. Robbiano and Zuily [21] consider smooth asymptotically flat nontrapping metrics 
in M" of the short range type and use a parametrix which is a Fourier integral operator with 
complex phase, relying considerably on Sjostrand's theory of the FBI transform. Hassell- 
Tao-Wunsch [11] have a more direct parametrix construction emulating the model of the 
constant coefficient fundamental solution, which applies to smooth asymptotically conic 
manifolds with short range scattering metrics, extended shortly afterward to long range 
scattering metrics. 

The dynamics for low frequencies are even more delicate, and for now there seem to be only 
two cases where anything at all can be said. The first is for sufficiently small perturbations of 
the flat metric, which is the case studied in [37j and here. The second is for time independent 
operators, with suitable spectral assumptions; for the Schrodinger equation this problem is 
considered in [17] (see also [2] and [22]), while for the wave it will be explored in another 
forthcoming paper. 

A key part of the global decay estimates are the local energy estimates, which measure 
the local averaged decay in the settings. In the simplest form (see e.g. [1], [12], [H], [16] . 
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[S], [20], [27], [30]), they are stated as 

when Dm = 0. Heuristically this is a reflection of the fact that waves move at speed 0(1) and 
thus spend a time 0{R) within a bounded spatial ball of radius R. These are the counterpart 
of the so called local smoothing estimates for the Schrodinger equation. See, e.g., [23], [39] , 
[1], [7], and [5]. A significant difference is that, in the case of the Schrodinger equation 
the speed is proportional to the frequency; therefore one also gains half a derivative in the 
estimates. 

The local energy estimates provide us with a convenient space to place the errors in our 
parametrix and also with a simpler setup in which to measure the decay of low frequency 
waves. In a nutshell, one of our main results asserts that 

Local energy estimates => Strichartz estimates. 

The most important part of the article is the outgoing parametrix construction, for which 
we are able to adapt the ideas in [37]. The parametrix construction in |37j is based on the 
use of a time dependent FBI transform. However it does not use Sjostrand's theory [24]. 
Instead, it takes advantage of the simpler approach introduced by the second author in [32j, 
[33]. 

For more information about phase space transforms, we refer to [9] and [^. One of the 
main starting points in the phase space analysis of pde's is Fefferman's article [H]. 

Simplified presentations of localized wave packet type parametrix constructions are now 
available in [15], [36]. These apply to evolutions of the form 

{Dt + a'"{t,x,D))u = 0, m(0)=mo 

on the unit time scale, for symbols a which satisfy a partial Sqq type condition 

\d:d^^a{t,x,0\<cap, \a\ + \P\>2. 

In the finite time analysis in [15], [36] the evolution is turned into a transport equation in 
the phase space modulo small errors. These parametrices are often useful in rescaled forms. 
However due to their finite time horizon, they cannot be directly applied to obtain optimal 
results for metrics which are not compactly supported perturbations of the identity. 

In the long time analysis in [37] a time dependent FBI transform is used instead. A 
second order term in an asymptotic expansion becomes nontrivial, and the equation turns 
into a degenerate parabolic evolution in the phase space. Bounds for this evolution are then 
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obtained using the maximum principle. Fortunately for us, the main step in the parametrix 
construction in [37] can be applied directly here for half- waves. See Theorem [291 

Even though our parametrix is very precise, there are still errors which need to be con- 
trolled and this is done using localized energy estimates. We prove such estimates in the case 
of small perturbations of the flat metric. For large perturbations nontrapping may fail, and 
thus the localized energy estimates may fail. A nontrapping assumption would help with 
the localized energy estimates at high frequencies, but not for the low frequencies. Here 
we avoid this problem by using the localized energy estimates as an assumption for large 
perturbations of the fiat metric. In the case of the Schrodinger equation, the local smoothing 
estimates for large perturbations were considered in [17]. See, also, [22]. In a follow-up paper 
we will consider the same issue in the case of the wave equation. 

Scaling plays an essential role in our analysis. Modulo rescaling and Littlewood-Paley 
theory all our analysis is reduced to waves which have fixed frequency of size 0(1). Since 
waves have a propagation speed of size 0(1), our study of outgoing waves can be largely 
localized to cones of the form ~ \t\}. Certainly the exact flow cannot have a precise 
localization of this type due to the uncertainty principle. To compensate for this we introduce 
an artificial damping term which produces rapid decay of waves which do not have the 
above localization. This allows us to restrict our attention to the above cone modulo rapidly 
decreasing errors. 

In the present article we consider global in time parametrices and Strichartz estimates 
for O^'^ metrics in R" which satisfy a weak asymptotic flatness assumption. Due to the 
global nature of the result it is convenient to consider scale invariant assumptions on the 
coefficients. We denote 

ylj = R X {2^ < \x\ < 2^'+^}, A<j = M x < 2^}. 

Following [37], we assume that 

(8) J^sup \x\'^\V'^a{t,x)\ + \x\\Va{t,x)\ + \a{t,x) - Mi+„| < e 

where Mi+„ is the (n + 1) x (n + 1) matrix diag(— 1, 1, . . . , 1) and, for the lower order terms. 




(10) sup |a;| |c(t, x)| < e. 

RxM" 
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In some special cases we will need to strengthen the last condition to 
(11) ^sup|a;|%(t,x)p < e. 

If e is small enough then ([8]) precludes the existence of trapped rays, while for arbitrary e it 
restricts the trapped rays to finitely many dyadic regions. 

Before we state our main results we need to introduce the function spaces for the localized 
energy estimates. We consider a dyadic partition of unity in frequency, 

oo 

1 = ^'^(^-)' 

fc=— oo 

and for each A; G Z we measure functions of frequency 2^ using the norm 

\\u\\xk = 2^/^||M||i2(^^_^) + sup |||x|"^n||L2(^^). 

j>-k 

To measure the regularity of solutions to the wave equation, we use the global norm 



|2 

fc=— oo 



All Schwartz functions -u G iS(M x M") have finite X'^ norm. This allows us to define the 
space as the completion of iS(M x M") with respect to the X* norm. Its structure is 
clarified by the next lemma: 

Lemma 1. [37] a) (s = 0) We have 

(12) sup |||2;|~^u||l2(^^) < \\u\\xo. 



b) If < s < then the following Hardy type inequality holds for all u G SiM. x R").- 
(13) |||x|~5-*^x||^2 < \\u\\xs. 

If ^ s < then we have the weaker bound 

oo 
j=-oo 

where the time dependent function ua^^ stands for the spatial averages of u in {\x\ < 2-'}. 

The proof of the lemma is similar to the special case s = | considered in [37j and is 
omitted. From the lemma we conclude that if s is as in case (a,b), then one can think of X 
as a space of distributions. On the other hand if s is as in case (c), then X has a BMO type 
structure, i.e. X is a space of distributions modulo time dependent constants. 



Controlling the constants is important, particularly when it comes to localizing paramet- 
rices in dyadic regions. This is why we introduce also a stronger norm which removes the 
BMO structure, namely 

iiii9 iiii9 Mil— 1 — cii9 ^ 

Iklljf. = Ikllx- + llkl u\\i2, 0<s<^— . 

This coincides with the X'^ norm for < s < To simplify the exposition we also set 
XO = XO. 

For the inhomogeneous term in the equation, on the other hand, we use the dual space 
= {X~^y with norm 

A;=— oo 

As X" is the completion of S{R x M"), for s > the space is dense in S'{R x M"). 

In addition, 

1 1 — Ti 

(15) \\u\\ys <\\\x\^-'u\\l2, < s < 
and 

(16) \H\yO <'^\\\x\^u\\L2(AJ)■ 

j 

Definition 2. We say that the operator P satisfies the H'^ localized energy estimates if for 
each initial data (mo,Mi) G H^^^ x and each inhomogeneous term f G L^H^ + Y'^, there 
exists a unique solution u to ^ with Vm G L°°H'^ fl X* which satisfies the bound 

(17) l|VM||^^^.nX= < ||Vm(0)||^. + \\f\\LiH^+Ys- 

In this context the lower order terms can be often treated as negligible perturbations: 
Lemma 3. a) Let b be as in (Q and 

Tl — 1 

(18) |s|<l, 1^1 
Then 

(19) \\bVu\\ys <e\\Vu\\xs. 

b) Let n > 3, c be as in f fTOj) and — 1 < s < 0. Then 

(20) \\cu\\ys <e\\Vu\\xs. 

c) Let n > 4, c be as in (fTTj) and s = —1, 0. Then 

(21) ||cM||ys < ellVwIlx"- 
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d) Let n = 3 and c be as in (|TT1) . Then 

(22) ||cM||yo < e||M||^i. 

The localized energy estimates hold under the assumption that the coefficients a"^ are a 
small perturbation of the Minkowski metric. 

Theorem 4. Assume that the coefficients a"^, 6" satisfy ([8]), iQ with an e which is suffi- 
ciently small. Assume also that c = 0. Then the operator P satisfies the localized energy 
estimates globally in time for s as in (ITSll . 

A general coefficient h and a coefficient c can be dealt with perturbatively but only in 
dimension > 3: 

Corollary 5. a) Let n > ?> and a°'^ , 6" and c as in 1^, IQ, ffTOl) with an e which is 
sufficiently small. Then the operator P satisfies the localized energy estimates globally in 
time for 

-1 < s <0. 

b) Let n > 4 and a""^ , 6" and c as in (fTTj) with an e which is sufficiently small. 

Then P satisfies the localized energy estimates globally in time for s = —1, 0. 

Once we have the local energy estimates, the next step is to construct an outgoing 
parametrix which has good time decay and suitable error bounds in the dual local energy 
spaces. The parametrix is constructed at ffist in the case of a small perturbation of the flat 
metric. This leads to our main scale invariant Strichartz estimate: 

Theorem 6. Assume that c = and the coefficients a"^, 6" satisfy ([8]), (Q with an e which 
is sufficiently small. Let {pi,Pi,qi) and {p2,P2,(l2) be two Strichartz pairs and s as in (ITSl) . 
Then the solution u to satisfies 

(23) ||VM|||i)^|Pi-.iPii,inx» < I|Vm(0)||^. + \\f\\^j^^^-,^-sLP'2L'''2+Ys- 

A zero order term c can also be added to P subject to the conditions in Corollary 

If e is large then any localized energy estimates require an additional nontrapping condi- 
tion. Even then the nontrapping can at most guarantee local in time bounds. However, we 
can still prove a conditional result: 

Theorem 7. a) Assume that c = and the coefficients a°'^ , 6" satisfy ([8]), (Q. Then for 
every Strichartz pair {p,p,q) and s as in ffTSl) . we have 

(24) \\\D^\'~^Vu\\lpli < ||VM||j^»nioo/fs + \\Pu\\y^. 
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In addition there is a parametrix K for P which satisfies 

(25) ||V/^/|||B,|p.-.^..i,.nx= + ll^/llx^+i + \\{PK - I)f\\ys < ||/|||^^|-.,-.^.^^,^ 

for any two Strichartz pairs (pi,Pi, qi) and {p2,P2, ^2)- ^ zero order term c can also he added 
to P subject to the conditions in Corollary {5^ 

b) Assume that in addition the operator P satisfies the H'^ localized energy estimates. Then 
the solution u to ^ satisfies the full Strichartz estimates in fl23|) . 

In applications one might be concerned that the condition ([8]) imposes the nontrivial 
restriction a(t, 0) = Mi_|_„. This is true, but it is needed only because we are allowing the 
derivatives of the coefficients to be singular at 0. Otherwise, such a restriction is unnecessary: 

Remark 8. Assume that the the condition ([8]) on the coefficients a"^ is modified for \x\ < 1 
to 

sup(|V^a(t,a;)| + |Va(t,a;)| + |a(t, x) - Mi+„|) < e, 

|x|<l 

and similarly for ([9]), ( fTOl) . and ( JTTI) . Assume also that for k > the definition of the space 
Xk is changed to 

IklUfc = \\u\\l^A<^o) +^^P\\\^\~^'^\\L^Ai)■ 
Then the results in Theorems^BS^ remain valid. Their proofs are essentially identical with 
only a few obvious changes. 

The paper is structured as follows. After introducing some notations in the next section 
and making a reduction to the case = — 1 in the third section, in Section H] we consider 
the paradifferential calculus associated to our problem. More precisely, we show that without 
any loss we are allowed to mollify the coefficients a"^ on a suitable x dependent scale. This 
allows us to reduce our analysis to problems which are frequency localized in dyadic regions. 
We also prove the bound. Lemma [3l for the lower order terms. 

Section [5] contains the proof of the localized energy estimates in Theorem HI The main 
step of the proof is carried out in a frequency localized context and involves a Morawetz type 
multiplier technique. 

After making a reduction to the half-wave operator in Section 6, we state our main result 
on the existence of frequency localized outgoing parametrices for half-wave equations, namely 
Proposition [15] in Section 7. Using this result we conclude the proof of Theorems [6117! 

The rest of the paper is devoted to the parametrix construction. This largely follows [37j. 

In Section [8] we introduce the pseudodifferential operators and the phase space transforms. 
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An important role is played by the conjugation of pdo's with respect to the phase space 
transform, for which we use some results from [M], [36]. In a first step, the parametrix is 
obtained in Section [9] in the case of evolutions governed by a pseudodifferential operator a}" 
whose symbol satisfies a suitable smallness condition uniformly in x; for this we are fortu- 
nately able to apply directly the result proved in [37]. This construction is then transferred 
in Section [10] to small perturbations of half-waves via conjugation with respect to the flat 
half-wave flow. Finally to arrive at the desired setup we need to insure that the parametrix 
is localized in outgoing propagation cones. This is done in the last section by means of 
choosing a suitable damping term in the equation. 

Acknowledgement: The authors thank the anonymous referee for a careful reading of the 
original version of this article and for the astute suggestions that have greatly enhanced the 
exposition. 

2. Notations 

We consider a smooth spatial Littlewood-Paley decomposition 

oo 
j=-oo 

We also set 

X<3 = Yxk- 

k<j 

Given e as in ([8]), we can find a sequence ej G so that 
(26) sup \x\'^\V'^a{t,x)\ + |a;||Va(t,a;)| + \a{t,x) - Mi+„| < ej 

A, 



and 



Without any restriction in generality, we can assume that is slowly varying, say 
(27) |lnej -lnej_i| < 2-1°. 

We also choose a function e in IR+ with the property that 

tj < t{s) < 2tj for 2^' < s < 2^+\ 



and so that 



This implies that 



|e'fs)| < 2-^s-ie(sl 



s 
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We also define ek{s) so that 

efc(s) ^ ej, s ^ 2\ j > -k 
efc(s) ^ e_fc, s < 2-K 

Note that 

efc(|x|)^e(2-^ + |a;|). 
We consider a frequency Littlewood-Paley decomposition 

oo 

1 = ^ S,{D,) 

j=-oo 

where 

supp C {2^-^ < lel < 2^+1}. 

We also use the related notations S^k, S^k, etc. 

We say that a function / is localized at frequency 2^ if / is supported in {2^^^^ < |^| < 
2'^'^^}. An operator K is localized at frequency 2^^ if for any / localized at frequency 2^^ its 
image Kf is frequency localized in {2'^~^° < |,^| < 2'^+^°}. 

3. A MINOR SIMPLIFICATION 

The aim of this section is to reduce the problem to the case when P has the form 

P = -Dl + 2Dia'^Dt + D.a^Wj + h'^D^ + c, 

and once this is accomplished, P will be taken to be of this form throughout the sequel. 
To arrange that = — 1 we multiply the operator P by — (a°*^)^^ which satisfies the same 
bounds as oP^. This modifies the other coefficients 

IP ^ -60(aOO)-i - (a°°)-i(Aa°°) + Dj{{a^')-')a^^, c ^ -c(a°°)-\ 

and it is easy to verify that the assumptions (E]), (E]), ( ITOi) . and ( fTTj) are left unchanged. 
To express the second term in the form above we note that 

This changes the coefficients 6° but still within the allowed limits. Arguing similarly and 
picking up only lower order errors within the limits, we may assume that a*-' = a^^. We also 
note that the coefficient c is not affected by these transformations. 

To conclude our simplification, we need to verify that our function spaces are not affected 
by multiplication by (a*"^)"^. 
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Lemma 9. Let a be as in ([8]) and s as in f|T8l) . Then 

(28) ||a/||y. < 

In addition, for all Strichartz pairs {p,p,q), we have 

(29) II«/II|D,|-P-=LP'L9' < \\f\\\D^\-P-^Lp'Li'- 

Proof. We write (128|) in the dual form 

\{af,u)\<\\f\\ys\\u\\x-s 

and take a simultaneous Littlewood-Paley decomposition of the three factors a, f and u. 
Nontrivial output is obtained when the two larger frequencies are comparable. Hence there 
are three cases to consider. The trivial one is when the a factor has the low frequency. For 
the remaining two cases it suffices to prove the off-diagonal decay 

\{S,aSkf,S,u)\ < 2(^-^)('=-^-)||5,/||xJ5,n||x,, 3 < K 

respectively 

\{S,aS,f,Sku)\ < 2^-^''^^'-^^\\S,f\\x'J\Sku\\x,, j < k 

for s as in ( JT8l) . This follows from the definition of the and Xk norms combined with the 
bound on Ska, 

\Ska{x)\<2-'\2-'' + \x\Y\ 
and an uncertainty principle bound for the low frequency factor on the dual spatial scale, 

\\SjU\\L2L^(^A^_^) < 2^^115-^.^11^^, \\Sjf\\L2L^(A^_^) < 2^mSjf\\x'y 

The details are straightforward and are left for the reader. 

We now prove (1291) . The time variable plays no role and is neglected in the sequel. We 
shall use the following variant of a Moser estimate: 

\\f9\\ws,p ^ ll/IU«i|klliy«.92 + II^I|l'-i||/IU=,-2 
with s > 0, 1 < p < oo, and 

- = \ = \ , g2, ^2 e (1, oo), gi, ri G (1, oo . 

p qi q2 ri r2 

See, e.g., [381 §2.1, Proposition 1.1, p. 105]. We first assume that s + p > and apply the 
above estimate to af. This yields 

(30) ||a/|| WP+!>,q' ^ ll'^l|oo||/||iyp+s,9' + ||c^|lvt'P+s,'-i ||/||r2- 
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The first term on the right is trivially bounded by the right side of fl29l) . For the second 
term on the right, we first pass from the Sobolev space to an appropriate Besov space, and 
then we use the following consequence of ([8]) 



^ \ |a.| \ ^ , m + / > 



\Sia{t,x)\ < 
Indeed, we see that 

k 

k m>—k 

Using that p = — for a sharp Strichartz pair, we have that q' < Thus, we may 
choose ri = Substituting this in the previous calculation, noting that p + s < 2 for p a 
part of a sharp Strichartz pair as above and s as in (ITSi) . and using the summability of {cm}, 
we have that 

n 

\\a\\wp+''--i < '^i = — : — • 

p + s 

If we now apply Sobolev embeddings to the second factor, we see that the second term in 
the right of fl30|) is also bounded by the right side of (|29|) . As we may use a dual argument 
if s + p < 0, this completes the proof. □ 



4. The paradifferential calculus 

In order to reduce the problem to a frequency localized context and to simplify the 
parametrix construction it is convenient to localize the coefficients in frequency. This is 
somewhat more complicated than usual because the frequency localization scale needs to 
depend on the spatial scale. 

It suffices to work with only the principal part of the operator P, which we denote by 

= -Dl + 2Dia'^Dt + D.a'Wj. 
Given a frequency scale k we define the regularized coefficients 

Kk-4 

Correspondingly we define the mollified operators 

14 



which are used on functions of frequency 1^ . Roughly speaking, their coefficients are fre- 
quency locahzed in the region 

lel <2'=(l + 2'^|x|)-i 
We also introduce a global mollified operator 

oo 

P = ^ P(k)Sk- 

k=—oo 

Due to fl26l) and to the fact that the e^-'s are slowly varying, it follows that the dyadic 
parts of the coefficients will satisfy the bounds 

2-2m-2l^^ 2m ^ |^| ^ 2m+\ m + l>0 

e_i \x\ < 2~'. 

This also allows us to obtain bounds on the coefficients of P{k), 

|a"(a!f Jx) - 6'^)\ < c^ek{\x\)2\^\\l + 2'=|x|)-l"l, |a| < 2 

(32) I 

< c„efc(|x|)2H'=(l + 2'\x\)-'--, \a\ > 2. 

The main result of this section shows that we can freely replace Pa by P in Theorems HI 
inilTl^a). It also shows that at frequency 2*^ the operators P and P(fc) are interchangeable. 

Proposition 10. Assume that the coefficients a satisfy ([H])- Then 

(33) ||(P - Pik))Siu\\x', < e||5/Vw|U„ 1/ - A;| < 2 



(31) 



l^^a'''(^,^)l < 



(34) mk),Sk]u\\x'^<e\\Vu\\x,. 
In addition, for s as in f ll8l) the following estimate holds: 

(35) ||(P„-P)^i||y. <e||VM||x- 

Proof. The proof is very similar to the analogous one in [HT]. We begin with fl5Sl) . and write 

Pa P Plow ~l~ Pfnid ~\~ Phigh 
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witli 

Plow 

Pmi. 



mid 



oo oo 

=— oo l<k—i k=—oo l<k—A 

oo k+A oo fe+4 



fc=— oo l=k—A 
oo 



k=—oo l=k—i 
oo 



P; 



A;=— oo i>fe+4 fe=— oo l>k+A 

Let us examine in detail the second term in each of the expressions above. The bounds for 
the remaining terms follow from similar arguments. 

For Plow, we notice that the output is at the same frequency 2^ as the input. Since the 
factor Dj contributes a factor of 2^, it suffices to show 



(36) I 

Kk-A 

Here, we shall use the bound 



iS<ix>k-2i)Sia^^v 



XL 



\S<iX>k-2iix)\ < 



2Al-Ak 

1, |a;|>2 



l<k-4. 



For Ixl ~ 2™, m > —k, we use this and fl3Tl) to see that 



— m— 1 



m— 4 



l<k—A /=— oo l=—m J— k-m 

'~ 2 

^ r\—m—k 



2m-2i 



For |x| < 2 the argument is easily modified to give 



Kk-A 



^ ^ — ky 



\x\ < 2" 



which, combined with the previous estimate, yields the desired bound 

For input frequency 2'^, Pmid permits output frequencies 2'^ for all h < k + 4. We take 
/ = k for simplicity of exposition, and consider separately low and high dimensions. 



^To be explicit with the order of operations, the Littlewood-Paley projectors take precedence, followed by 
multiplication and differentiation using right associativity, and finally addition. Thus, for example, the first 
term in PiowU is understood to be 



oo 

E S iS<iX>k-2i){Sia'')){D,{Sku)} 



k— — c 



Kk-i 
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In low dimension n = 2,3 the bound for Pmid follows from the off-diagonal decay 
\\ShDi{Ska''DtSku)\\x'^^ < e2^('^-'=) || A^fc^llx,, h < k + 2, 

or more simply, 

(37) \\S,{Ska%)\\x'^^ < e2'^^^-'h~%\\x,. 
Writing 

m>—k 

it is sufficient to show that 

\\Sh{XmSka'%)\\x'^ < e^2'^^^-'h~''\\v\\x„ m > -k, 
\\SHix<-kSka''v)\\x'^ < e^,2'^^''-'h-'\\v\\x,. 

For the former, we apply ( |3T1) and see that it suffices to show 

(38) \\Sh{Xn.v)\\x', < 2'^('^-'=)2^+'?'||t;|U2(|,|^2'"), m + k > 0. 
By interpolating the estimates 

\\Sh{XmV)\\L^ <\\XmV\\Lh \\xS h{XmV)\\ < 2"'\\XmV\\ , TU + h > 0, 

we obtain 

(39) \\ShiXmV)\\x',^<'2^\\XmV\\L^- 

Recalling that we are in the case when h < k — 2, this yields (1551) when m + h > 0. 
For m + /i < we have improved bounds 

n{7n-\-h) I n{m-\-h) 

\\ShiXmV)\\L2 < 2^^||x™^^||l2, \\xShiXmV)\\L2 < 2-^2^^ \\XmV\\2, TU + h < 0, 

which upon interpolation yields 

l|5/.(Xn^^^)llx; <2t2^(™+^)||x™t;||L^, 

and implies (!38|) . The bound when Xm is replaced by x<~k is identical to the m + h < 
argument above. 

In high dimension n > 4 the bound flTTI) is replaced by 

(40) \\ShiSka'%)\\x', < e2'^'^-''>2-'^\\v\\x,. 

whose proof is similar. The only difference is that now the worst case is m = —h, whereas 
in low dimension the worst case is m = —k {n = 2) respectively —k < m < —h (n = 3). 
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It remains to consider Phigh whose output is at frequencies 2' with I > k where 2^ is the 
input frequency. In low dimension n = 2, 3 it suffices to show that 

WSiD^iSia'^^DtSkU^xi < e2^('=-')|| A^fc^Hx,, l>k + 4 

or 

(41) WSiiSia'^'SkV^xi < e2'^^'-^^2-'\\Skv\\x,, l>k + A 

which follows by duality from (!37|) . 

In high dimension n > 4 instead of fj^Tl) we have 

(42) \\Sia''Skv\\xi<e2''-'2-'\\v\\x, 

which is still sufficient except for the endpoint s = 1. At the endpoint we are left with no 
off-diagonal decay. To compensate for that we need a stronger version of fH2]) . namely 

(43) WSyk+ADy'^SkvllYO < e.k2'\\v\\x,. 
By f ITBl) it suffices to show that 

oo 

\\\x\^Syk+iDla'^Skv\\L2(^A,) < e-fc2^||^^IUfe- 

j=-co 

But this follows from the bound 

\S^k^,Dy'\<e,{\x\)\xr', 
and in the case that j < — /c, the following consequence of Bernstein's inequality 

\\X<-kSkV\\L2L^ <2'^''\\SkV\\xk- 

The constant e_fc is obtained since the worst case is when j = —k, with exponential decay 
away from it. 

The estimate (133!) follows from arguments similar to those used for Plow To show it 
would suffice to show 

and the equivalent statement with j = 0, which follows directly from (1321) with a = 1. □ 

In a similar manner we prove the bounds of Lemma [3] which show that in high dimension 
we can completely dispense with lower order terms. 

Proof of Lemma O This is again quite similar to the related result from [37J . 
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From iQ we may obtain the following bounds on the frequency localized pieces of the 
coefficients 

(44) \Skb{x)\ < 2-hk{\x\){2-' + \x\)-\ \S^kb{x)\ < ek{\x\){2-'' + l^r^ 

To prove (fT9i) . we first expand 

j j |fc-i|<4 j fc>j+4 

The easiest case is the first term: the low-high interactions. Here the output is at the 
same frequency range as the input. Thus, it would suffice to show 

\\{S<j^ib)v\\x'^ < 4v\\xy 

As the e(|x|) provides summability, this follows directly from (jS]). 

The third case is the high-low interactions. Here, for input at frequency 2^ , the output 
is at frequency 2^ with k > j + A. We, thus, measure the output in X^. In low dimension 
n = 2, 3 it suffices to show that 

\\SkbSjv\\xi < e2^(-'-'=)||5,t;|U,, J + 4 < fc. 

This, however, is just a reformulation of fj4T]) since b has exactly the same regularity as Da. 
In high dimension n > 4 we have the similar relation 

\\SkbS,v\\x'^ < e2^'-^|| Vllx,, j+4<k 
which covers all cases but s = 1. For s = 1 we replace this with (H3i) with Db instead of 

The remaining case, the high-high interactions, is dual to the previous case. This finishes 
the proof of dUl). 

Finally, (EOl), ([21]), and ([22]) follow directly from the embeddings ([12]), ([13]) and their duals 
([ED- □ 

5. Localized energy estimates 

Here we prove Theorem |H We can assume that P has the form in Section [3] with c = 0. 
Also due to ( |T9]) we can take b = 0. 

The theorem is proved via a positive commutator method. Let (am)m6Z be a positive 
slowly varying sequence with = 1- Correspondingly we define the space ^ with 

norm 

j>-k 
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and the dual space 

j>-k 

The key step in the proof of Theorem H] is the following frequency localized estimate: 
Proposition 11. Assume that e is sufficiently small. Then the hound 

(45) llVKllL-L^nXfc,, < ||Vn(0)||L2 + \\P{k)u\\L^L^+x'^^^ 

holds for all functions u G L°°L^ fl Xk^a localized at frequency 2^ , uniformly with respect to 
all slowly varying sequences (am) (as defined in (1271) ) with 

oo 

(46) 5]; «™ = 1. 

m=—k 

Proof. By rescaling, the problem reduces to the case when k = 0. We may without loss 
increase the sequence {am) so that it remains slowly varying and 

This is accomplished by redefining 

e 

Since e <C 1, we can fix another small parameter e ^ 5 ^ 1 so that 

Associating to (a^) a function a{s) = ao{s) whose definition is analogous to that of €q{s) 
in Section [2], we have, from the last property, 

^o('S) ^ ea(s) ^ 6a{6s). 

The proof has three ingredients, the first of which is the classical energy estimate. Since 
e is small it follows that the dt vector field is time-like, and the corresponding energy 

E,iu) = ^\\DM\' + l{a%D,u,D,u) 

is positive definite. Here and throughout (■, ■) is the L^(M") inner product. The time 
derivative of this energy is 

j^Eoiu) = 53(P(o)M, Au) + ((9,afo))Aw, Aw) + ^((^tajy^.M, A^)- 
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The second component of the proof is a Morawetz-type commutator estimate. Let Q{x, Dj.) 
be a spatially self-adjoint operator. On time slices we obtain 

(48) ^|-23?(A^i, Qu) + 25R(ajjD,M, Qu)'^ = -2%P^^)U, Qu) + {i[D,al^D,, Q]u, u) 

+ 23f?(2[afo) A, Q]u, Dtu) - 2n{{d,a%^)DtU, Qu) + 2K((ata;°)) A.m, Qu). 

The point here is that we seek to choose Q so that the commutator [Dia^^^^Dj, Q] is positive 
on the characteristic set of the operator P(o). 

Finally, to account for the elliptic region, i.e. away from the characteristic set of the 
operator P(o), we use a Lagrangian term. Precisely, for a real- valued, time-independent, 
scalar function ip{x), we compute 

^^{{-Dt + 2a°^^^Dj)u,^pu) = 3fJ(P(o)U,^/'u) - ^{a^Q^DjU,^/jDiu) + ^{a^^^^DjU, {di^lj)u) 

+ 2$5((ata°j)D,M, iju) - 2!5((a,a°^)) An, i^u) 
+ (An, ^IjDtu) - 23?(a°j A-M, ^IjDtu). 
We consider two additional small parameters 5q and 5i so that 

and define the modified energy 

E{u) = Eo{u) - 5o3f?((-A + all^Dj)u, Qu) - (5i3?((-A + 2a°^'^ A)«> 
Combining the last three relations we obtain 

-E{u) + -j{t[D,a'^^Dj,Q]u,u) + 6i{Dtu,tpDtu) < 3(P(o)U, (A + SqQ + iSiij)u) 
+ {\Va^o)\^u,Vu) + 6o\{i[al^)D,,Q]u,Dtu)\ + 6o{\Va^o)\\Vu\, \Qu\) 
+ SM^xu\,^\Vu\) + 5i(|a||V,M|, \Vij\\u\) + (5i(|Va(o)||VM|, \i;\\u\). 
We choose Q as in [27]. For convenience its properties are summarized in the following 

Lemma 12. There exists an operator Q of the form 

Q{x, A) = 6{Dx(j){6\x\) + (f){6\x\)xD) 

where has the properties 

(i) 0(s) ^ (1 + s)-^ fors>0 and \d^(f){s)\ < (1 + s)~''~^ for k < 4, 

(li) (f){s) + s(f)'{s) ^ (1 + s)"^a(s) for s > 0, 

(Hi) (p{\x\) is localized at frequency <^ 1, 
and which satisfies the bounds 

\\Qu\\l2 < ||m||l2 
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{i[Dia,l.Dj, Q]u, u) dt > 5\\ 



for all functions u localized at frequency 1. 
The function '?/'(|a;|) is chosen so that 

1 + s 

We first note that the above properties of Q and ip insure that E is positive definite; specif- 
ically 

E{u) ^ \\Vu\\l. 

for all functions u at frequency 1. Moreover, upon integration in t, we can estimate 

{D,u,ijDtu)dt>\\DM\\^. 
and thus, the integral of the left side of fH9|) is bounded below by 

snvE{u){t)-E{um + 5\\VM\\,^+5Mn\\\,^- 
We now examine the right side of fj49l) after integration in t. Using (!32|) . we have 

(|Va(o)|Vn, Vu) + \{i[af^.^D,,Q]u, Dtu)\ + (|Va(o)||Vu|, \Qu\ + \ij\\u\) dt < e||VM||^^_^. 
Similarly, by our choice of ip, "we may find a constant M > so that 



j C6M\VM,i^\'^^\) + C6i{\a\\VMA^i^\u)dt<j\\dM\^^^^^ 



where C is the implicit constant in fH9l) . 

Using these bounds to estimate the right side of and using Cauchy-Schwarz, we obtain 

provided, say, 66o > 2M6i. This concludes the proof of Proposition [TTl □ 

We conclude now the proof of Theorem HI Let be another slowly varying sequence 
with 

m 

Applying Proposition [11] with am replaced by am + Pm we obtain the bound 

||VM||L-L2nXfe,,+^ < ||Vm(0)||l2 + II P(fc)M 11^1^2+^/^^^^ 

for all u localized at frequency 2^. This implies the weaker estimate 

||Vn||i^i2nXfe,, < ||Vu(0)||l2 + \\P{k)u\\L^L2+x'^_^- 
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Since any sequence is dominated by a slowly varying sequence, we can drop the assump- 
tion that a and (3 are slowly varying. Then we maximize the left hand side with respect to 
a G and minimize the right hand side with respect to /3 G This yields 

(50) ||VM||Looi2nXfc < ||VM(0)||i2 + \\P(k)U\\L^L2+x'^. 

For an arbitrary function u G X^, we apply this bound to SkU. We have 

P{k)SkU = SkPu + [P(^k),Sk]u + Sk{P{k) - P)u. 
The last two terms are frequency localized and can be estimated by (1331) and ([3i 
mk),Sk]u + Sk{Pi^k)-P)u\\x',<t J2 W'^SiuWx,. 

\k-l\<2 



Then after summation we obtain 



k 

'^umHs + J2[^"''\\SkPu\\i.,.^x', 



< 



+ 2''''\\[P^k),Sk]u + SkiP^k)-P)u\ 



< l|V«(0)f„. + IIP^II^,^.^^. + e|| Vm||^. 



< II V^i(0)||5,. + ||Pa^i||i,^.+^. + e||Vn||^. ( by m) 

< \\VuiO)\\% + \\Pu\\l,^^^y, + e\\\/u\\ls ( by Lemma[3]). 
For small e we can neglect the last right hand side term to obtain 

(51) \\Vu\\l^^,^^, < ||Vm(0)||^. + \\Pu\\l,^,^y, 

which holds in any time interval containing 0. 

Reverting the transformation in Sectional we see that without any restriction in generality 
we can write P in its self-adjoint divergence form. Assuming that 6 = 0, we may then use a 
duality argument to show that for any / G L^H'^ fl Y^, there is a f solving 

Pv = f, v{0) = Vo, vt{0)=vi 

with 

l|Vf ILoo^^nx^ ^ l|Vt;(0)||^. + 11/11 
Due to f|T9|) this extends perturbatively to the case of nonzero b. 

By (ISTl) . this solution is unique, and the proof of Theorem H] is concluded. 
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6. The half wave decomposition 

In this section we reduce the study of the wave equation ([7]) to the study of two half-wave 
equations. We first factor the principal symbol as 

-r^ + 2aOVe, + a'^^.^j = -(r + a+(t, x, 0)(r + a-{t, x, 0) 

where are 1-homogeneous in ^ satisfying the symmetry property 

a~{t,x,^) = —a^(t,x,—^) 

and are chosen so that > . The symbols can be written down explicitly as 



In the sequel, we shall, however, only need the properties listed above. This will permit us, 
in Section [H] and beyond, to free up the a, b, c notation. There the focus will only be on the 
half- wave operators and the symbols a^. The notations a, b, c will no longer be reserved for 
the coefficients of P but will be used for abstract terms which play the analogous roles. 

Mollifying the symbols with respect to x as in Section H] we obtain the symbols 
a^f^^{t,x,^) which we use at frequency 2^^. We note that are not the symbols obtained 
from the factorization of the principal symbol of P{k)', also one cannot define them in this way 
since algebraic operations (such as square roots) do not preserve the frequency localization. 

We also denote 



and let l(k)it, x,C,) be the corresponding regularizations. We note that l(k)(t, x,^) is obtained 
by regularizing /(t,x,^) and not by algebraically combining the symbols a^^(t, x,^). 

We are interested in operator properties matching the above algebraic properties. We 
work at frequency 1, but by rescaling the results extend to all dyadic frequencies. 

Proposition 13. Define the error operators 



(52) i?+ = P(o) + (A + A(o))(A + Af), R- = P(o) + (A + 4o))(A + A(o)). 




^(^,a;,0 = (a"^(^,a;,0 - a {t,x,^)) 



-1 



Then for all functions u and f localized at frequency 1, we have 



(53) 








(54) 



(x)(L(o)(A 



+ 



A7j-I)f\\Li< 




(0) 



(55) 



ll[^{0),^(0)]^^||x;^ 



< llVnllxo- 
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Proof. We write in the form 

R+ = - idtA+^{t, X, D) - 22(9,aJ° )A - 

- (A-A+)(o)(t, X, D) + (t, a;, a;, D). 

The first three terms are easily estimated by (!32|) . Consider the remaining two terms. For 
1^1 ^ 1 the symbols a^(t,x, ^) are smooth and homogeneous in ^. Expanding them into 
spherical harmonics we can assume without any restriction in generality that both have 
the form 

a^{t,x,i)=a^{t,x)h^{i) 
where a^(t,x) satisfy bounds similar to the bounds for a*-^, namely 

(56) \a^{t,x) - a^l + \x\\Va^{t,x)\ + \x\^\V^a^{t,x)\ < e{\x\). 

Then the last two terms in R'^ have the form 

a^o)(t,x)/i"(L')aJ,)(t,x)/i+(L') - (a"a+)(o)(t, = 

(a^o)(t,a;)aJ))(t,x) - {a~a+)(^o)it,x))h'{D)h+{D) + a^Q)(t, a+)(t, a;)]/i+(L'). 

The operators /i^ are bounded in Xq on functions of frequency 1. The commutator estimate 

(57) [h-{D),al,^{t,x)]:Xo^X;, 
on frequency 1 functions follows due to the bound 

|Va(Q)(t,x)| < eo{\x\){x)~\ 
Hence the estimate ( l53l) is proved if we can show that 

Lemma 14. Let a"^ be functions satisfying fl56l) . Then 

(58) |a(Q)(t,x)aj3)(t,x) - (a"a+)(o)(t,x)| < eo(|a;|)(a;)"\ 

Proof. Without any restriction in generality we can assume that = 0. As in the case of 
the coefficients the regularized functions have size 

|afo)(t,a;)| < eo(|x|). 

We separate the contributions coming from small x and from large x. The contribution from 
small X decays rapidly at infinity, 

\{x<oa^)(o){t,x)\ < eo(|a;|)(a;)"^, 

and the corresponding part in fl58|) will satisfy a similar bound. Hence without any restriction 
in generality we assume that are both supported in A>q. This allows us to replace (!56l) 
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with a better bound 

|a±(t,x)| + {x)\Va^{t,x)\ + {xf\V'a^{t,x)\ < eo{\x\). 
Using the analogues of (!3T|) . this allows us to estimate the differences 

Wto)it^^) - a^{t,x)\ < eo(|x|)(x)~^ 
and similarly for their product a^a~ . The conclusion of the lemma follows. □ 
The proof of ( l54l) is virtually identical, the roles of a'^itjX,^) are played by l(t,x,C,) and 

For fl55|l we expand /(t,x,^) in spherical harmonics and reduce the problem to the case 
when 

i{t,x,o = Kt,x)h{0 

with l(t,x) satisfying (156|) . Then the proof of (!55l) reduces to commutator estimates similar 
to (EZD. □ 

7. Parametrices and Strichartz estimates 

Here we reduce the proof of Theorem [6] to the construction of a suitable parametrix for 
Dt + ^foy Our main result concerning parametrices is 

Proposition 15. Assume that e is sufficiently small. Then there are parametrices for 
Dt + which are localized at frequency 1 and have the following properties: 
(i) L"^ bound: 

(59) \\K^it,s)\\Li^Ll<l, 
(a) Error estimate: 



(60) 



11(1 + \x\f{Dt + Af^^)K^{t, .)1U.^^. < (1 + |t - ^1)-^, t ^ s, 
[1 + \x\fDt{Dt + Af,^)K^{t, s)hi^Ll <{l + \t- 5|)-^, t + s, 



(Hi) Jump condition: K^(s+0, s) and K^(s—0, s) are type pseudodifferential operators 
satisfying 

{K^{s + 0,s)- K^{s - 0, s))So = So, 

(iv) Outgoing parametrix: 

(61) ||l||,|<2-io|i_,|}iro±(t,s)|U|^i| < (1 + |t- s|)-^ 

(v) Pointwise decay: 

(62) \\K^{t, s)hi^L^ < (1 + It - ^1)^"^- 
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Here is defined by 

/oo 
K^{t,s)f{s)ds. 
-oo 

We leave the proof of tliis result for later sections, and we show that it implies Theo- 
rems [6|71 As an intermediate step we have the following localized Strichartz estimates for 
the parametrix: 

Proposition 16. The parametrix given by Proposition [731 has the following properties: 
(i) (regularity) For any Strichartz pairs {pi,qi) respectively (^25^2) with qi < q2 we have 

(63) l|i^o^/llL.iL.inx., < 11/11^/,^,^ 
(a) (error estimate) For any Strichartz pair {p,q) we have 

(64) ||[(A + 4))i^o^-l]/||x^, < 

In both (!63l) and ( !64l) the function f is assumed to be localized at frequency 1 . 



The proof is identical to the proof of the similar result in [STJ Proposition 12] and is 
omitted. The proof of fl63l) follows that of the Strichartz estimates in the constant coefficient 
case as it consists of interpolating between fl59|) and fl62l) . using a TT* argument, and applying 
the Hardy-Littlewood-Sobolev inequality. The error estimate (16^ follows somewhat directly 
from (EQl). 

We can use the half-wave parametrices to construct a full wave parametrix. Precisely we 
have 

Proposition 17. Assume that e is sufficiently small. Then there is a parametrix Kq for P(o) 
which has the following properties: 

(i) (regularity) For any Strichartz pairs {pi,qi) respectively (p2, Q'2) with qi < q2, we have 

(65) ||V7ro/|U.iL.inXo < ll/IL/,^,^ 

(ii) (error estimate) For any Strichartz pair {p,q) we have 

(66) \mo)Ko-l)f\\xi,<\\f\\L.'L.'- 

In all of the above the function f is assumed to be localized at frequency 1 . 

Proof. Our first approximation for Kq is the operator Kqq defined by 
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The operator L(o) is bounded in both L^^L"^ and Xq] therefore from (l63l) we obtain part of 
(!65|) . namely 

We can also bound DfKoof in Xq. We have 

D.Koof = [A, ^(0)] W - ^o") + ko)Dt{K - Ko), 

and the first commutator is bounded in Xq. For the second term we use the Xq bound for 
L(o) and write 

Now we use ( l64l) and the embedding Xq C Xq for the first two terms and the Xq boundedness 
of and A^j. Summing up we have proved that 

We still have to estimate DtKoof in W^L'^^^ but we postpone this for later. 
Next we estimate the error 

The kernel KoQ^s^t) of Kqq is smooth in s, t away from the diagonal. However, we need to 
compute its singularity on the diagonal. Due to the property (iii) in Proposition [15] we see 
that the jump of i^oo on the diagonal vanishes, namely 

[K^oit, t)] := Koo{t + 0,t)- Koo{t - 0, t) = 0. 

However, the jump of the t derivative of KoQ{t, s) on the diagonal is nontrivial. Precisely, we 
have 

[DtKoo{t,t)] = L(^o)[Dt{K+ - K-Kt,t)] 

= L(o)(-A+)[Ko+(t,t)] + A-)[i^o-(t,t)]) 

+ L(o)([(A + Al^^)K^it,t)] - [(A + Al,^)K^it,t)]) 

= koMio) - ^J)) + ^(0)([(A + - [(A + A(o))i^o"(t,t)]). 

By ( IMI) the first term is close to the identity, while the second can be estimated by ( l60l) . 
For / localized at frequency 1 we obtain 

(67) 11(1 + |x|)([Ai^oo(t,t)] -/)/IL.^^2 < < ll/IL.^^.^ 

Next we compute 

(P(o)iroo - 1)/ = i?o/ + ([A^oo(t,t)] - /)/ 
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where the first term represents the off-diagonal contribution and the last term represents the 
contribution due to the jump of DtKQQ{t, s) on the diagonal. 

We use the factorization fl52l) for P(o) to compute the kernel of Rq, 

Ro{t, s) = L(o)P(o) W - Ko)it, s) - [L(o), P(o)] W " Ko){t, s) 

= - L(o) ((A + A(o))(A + A+^)K^it, s) - (A + 4o))(A + ^(o))^o"(t, s)) 

+ L(o) {R'-K+it, s) - R-K^{t, s)) - [L(o), P(o)] W " K,){t, s). 

For the expression on the first line we use the Xq boundedness of L(o) and A^^y together 
with the error estimates in fl60l) . For the R"^ terms we use f l53l) together with fl63l) and fl60|) : 
the latter is needed to bound the time derivative DtK(t,s). Finally, for the last terms we 
use (l55l) . Summing up, we obtain 

(68) \\Rof\\x!,<\\f\\L.'L.'- 

This is an acceptable error. 
The expression 

f, = -i[DtKoo{t,t)]-I)f, 

however, is not an acceptable error because it does not yield to a similar bound of its Xq 
norm. It has better decay at infinity; therefore we can account for it by setting 

Kof = Koof + Koifi 

where 

1 r°° 

KoMt) = -- / e-l*-^lA(s)rfs, 

which solves 

d^KoJ, = fi + KoJ,. 
Using the bound (1671) for /i, it is easy to see that iCoi/i satisfies 

II V-ft'oi/lllLPiL'JinX,, < ll/lliP'2i9^- 

On the other hand the /i component of the error is replaced by 

h = P(o)i^oi/i - /i = (2Aafo) A + Aag)^, + 1)^^01/1, 
which we can estimate by 

ll/2||x^ < llA^oi/ilU^ + ll^oi/ilU^ < 1 1 (a;) /i 1 1 ^2 
and then apply (!67|) . 
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The last step of the argument is to prove the L^^L'^^ bound for DiKq/. We will show that 
for u at frequency 1 we have 

from which the desired bound follows after an application of (!66|) . This would follow from 

II A^IIlpiLW ^ ||m||lpiL9i + \\P(0)U\\(L^+LPi)Ln 

or equivalently, 

llAwllLfiLn < \\u\\LPiLn + Wgih^Ln + ||5'2 1| Lf iL« , P{0)U = gi + 02- 

The above is easily reduced to the case (^i = by substituting 



u:=u — v, ^ ~ ~2 ' ^ 9i{s)ds 



oo 



smce 

||'^||(LinL°°)L9i + II A^||(LinL°=)L9i ^ ||5'i||lil9i- 

We are left with proving 

llAt.ll LPiLn ^ ||lt||LPiL'Ji + II^(0)M|| 
which follows from the interpolation inequality 

l|AM||ipiLw < ||m|UpiLw||A^m||lpiLW. 

□ 

Proposition [T7] is useful only if e is small. However, a similar result holds even if e is not 
small: 

Proposition 18. Assume that the coefficients a^^ satisfy (IHl)- Then there is a parametrix 
Kq for P(o) localized at frequency 1 and which satisfies 

(i) (regularity) For any Strichartz pairs {pi,qi) respectively (p2,i?2) with qi < q2, we have 



(69) \\VKof\\rv. ra,r.^ < 

(a) (error estimate) For any Strichartz pair {p,q), we have 

(70) ||[P(o)iro-l]/|U^ < 
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The proof is identical to the proof of the similar result in [371 Proposition 15] and is 
omitted. The idea is that the smallness condition is violated only on finitely many dyadic 
spatial regions. In [37] it is argued that a fixed dyadic spatial region can be partitioned into 
finitely many cubes on which the smallness holds with respect to a different coordinate frame. 
The local parametrices are then assembled together using a partition of unity. Alternatively, 
in a fixed dyadic region the problem of constructing a parametrix as above can be localized 
to a similar time scale and then rescaled into a local problem. 

Proof of Theorems 0, In what follows we work in a time interval [T^, possibly infinite. 
By (133|) we can replace the operator P(o) by P in Propositions [TTl [HI Rescaling this result we 
obtain similar parametrices Kj at any dyadic frequency 2^ . We first assemble these dyadic 
parametrices and set 

oo 
j=-oo 

The properties of K are summarized in the next lemma. 

Lemma 19. The parametrix K for Pa has the following properties: 

(i) (regularity) For any Strichartz pairs {pi,Pi,qi) respectively (p2,P2,?2) with qi < q2, we 
have 

(71) ||Vir/|||o^|Pi-.iPii,inx= ^ \\f\\\D^\-P2-sLp'2L''2- 
(a) (error estimate) For any Strichartz pair {p,p,q), we have 

(72) mK-I)f\\ 

y ^ 11/11 |D|-P-=LP'L9'- 

Part (i) follows directly from the Littlewood-Paley theorjl^. Similarly we get part (ii) but 
with P instead of Pa, since we can write 

PK-I = J2iP - PU))KA + {Pu)Kj - I)S, 

jez 



However, by ( ]35ll we can freely interchange Pa and P. Since (jlQll allows us to further pass 
from Pa to P with c = 0, this establishes the bounds for the first and last terms in the left 
side of ([25]). 

A second step is to use duality to establish an — > L^L^ bound. This establishes ([2] 
i.e. the first part of Theorem [7| (a). 



■^As mentioned in the Introduction, the Littlewood-Paley theory with respect to the spatial variables 
cannot be used in dimension n = 2 for the L^L°°, respectively the L'^/^L^ norms. Here we instead only 
obtain estimates in appropriate P Besov spaces. 
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Lemma 20. // there is a parametrix K for Pa as in Lemma [7PI and {p,p, q) is a Strichartz 
pair, then 



(73) 



\'^'^\\\D^\P-'LPLi ^ II Vull^oo^fsnX" + ll-Pfl'^lll'^- 



Proof. Without any restriction in generality we assume that and are finite but prove 
the bound with constants which are independent of T"*" and T~ . For g"" G {Dxl^'^'L^' L'^' we 
use integration by parts 



/•r+ /•T+ pT+ 

/ {Vu,g)dt= / {Vu,PaKg)dt- 
It- Jt- Jt- 



{Vu, [PaK-l]g) dt 



-{Pau, V-Kg)- {Vu, [PaK - l\g) - 2{d,u, (Va*°) ■ dtKg) 



dt 



- {dju, (Va^^) ■ d,Kg) + 2{diU, idta''')V ■ Kg) - 2{dtu, (a,a'°)V ■ Kg) 

+ {Vu, dtKg)\li,- + {dtu, V ■ Kg) - 2{d^d,u, V ■ Kg) |?! 

- {dtU,dtKg'')\'^l + 2{a^'d,u,dtKg'')\^l + {a^^d,u,d,Kg')\^l. 
Then by ( ITT]) and (|72]) we obtain 

J {Vu,g)dt <\\g\\\D\=-PLP'Li'(\\'^'^\\L--H=nX'^ + \\PaU\\Y'y 

Here we have also used f|T9|) with b replaced by Va, which according to ([8]) satisfies ([9]). The 
conclusion follows. 



□ 



Next we prove that the conclusion of Lemma [T^ is also valid for qi > q2: 

Lemma 21. The parametrix K in LemmalWi also satisfies (I7T1) when qi > q2- 

Proof. We repeat the computation in the previous lemma with 

u = Kf, g"" e\D\'-P'LP'^L'''K 
All the terms are estimated in the same way except for 

{Pau, V-Kg)dt= / {{PaK -I)f,V- Kg) dt + (/, V ■ Kg) dt 



T- JT- 

for which we use f lTTl) and fl72|) to estimate 

[ {PaU,V-Kg)dt <\\iPaK-I)f\\ys\\VKg\\ 
Jt- 



X- 



|_D|-s-P2L''2L92 



l|Vir^7|| 



< 



\D\-'^-P2LP'2L''2 II IDl^-PiL^iL"! ' 
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Then as in the previous lemma we obtain 

r.T+ 



{Vu,g)dt 



T- 



~ \\9\\\D\^-PiL^'iL< \ \\^^\\L'^H''r\X'' + 11^ ll|D|-''-P2Lf2L'J2 
~ Il5'll|£i|s-Pi7;^pij;^9i II/II|j:)|-s-P2L''2L''2 



which concludes the proof. □ 

The bound (!7T|) on Kf allows us to estimate HK/Hx^+i- However, if s + 1 > then in 
order to conclude the proof of (1251) . i.e. the remainder of Theorem [7] (a), we need to have a 
bound for the stronger norm This is achieved in the next lemma. 

Lemma 22. There is a parametrix K for Pa which satisfies: 

(i) (regularity) For an^Strichartz pairs {pi,Pi,qi) respectively {p2,P2,<l2), we have 

(74) \\'^Kf\\iD.in-sL..Lnnxs + ||A7b^^+i < 

(a) (error estimate) For any Strichartz pair {p,p,q), we have 

(75) ||(P„^-/)/||y. <||/||,^,_,_.^,,^,,. 



Proof. Let K be as in Lemma [I9l If we think of Kf as the sum of its dyadic pieces which 
are measured in X^, then for s + 1 > we fail to obtain a X'"^^^ bound for Kf due to 
the accumulation near the origin of the contributions below the uncertainty principle scale 
{|^l ~ 1^1^^}- To remedy this we attempt to remove these contributions. 
We consider a Schwartz function with 

0(0) = 1, supp0c{|e|G [1/2,2]} 

and set (f)k{x) = (j){2^x). In a first approximation we replace the parametrix K with (1 — T)fC, 
with T defined by 

oo 

Tu= TkSkU, TkU = u{t, O)0fc. 

fc=— oo 

This substitution improves the left hand side of ([71]). We shall show that 

Tl — 1 T? -I- 1 

(76) 11(1 - T)nb.,. < -^^' + ^<^ 

(77) WVTKfW 

For functions u localized at frequency 2^^, we have the fixed time pointwise bound 

Ht,o)\<^"^'\Mxo 



^We are again largely ignoring the L'^L°°, respectively L^/'^L^, estimates in n = 2. 
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which imphes that 

\\Tku\\x^ < \\u\\x^- 

Here, X° is the spatial part of the Xk norm, i.e., Xk = L'^X^. Hence we easily obtain 

The X^ bound of (1761) follows immediately, and in order to obtain the X** bound of (177|) . 
we then apply ( 1711) . The U'^L'^^ estimate uses a similar argument involving a Bernstein 
estimate, Littlewood-Paley estimates, and the bound flTTj) . 

To prove the part of the bound (!76l) . we take advantage of the fact that ((1 — 
Tk)Sku){t, 0) = to obtain the better bound 



_-|_n 

sup|||a;| ^{\x\+2~'')^+-{l-Tk)Sku\\L2(A,) <\\Sku\\x„ 
j 

which after summation yields 

1,1 n ^ n n 71—1 U+l 

\\\x\ 2(1 -T)M||i2 < + 
Consider now the error estimate for {1 — T)K. We claim that 

(78) \\PTKf - TfWvs < \\{Pk - l)f\\ys + WVKfWxs. 

It is easily seen that Tk is bounded in X^; therefore T is bounded in Y^. It remains to show 
that 

\\PTu-TPu\\ys < \\Vu\\xs 

which reduces to 

\\P(k)TkSkU - TkP(k)Sku\\x'^ < \\VSku\\xt,- 
After rescaling to A; = this is straightforward. What is important is that the second order 
time derivatives cancel. All the remaining terms can be estimated separately. 

It remains to consider separately the outstanding error estimate for TKf. This cannot 
be placed in because it does not have enough time integrability. Hence we need to add a 
correction to the parametrix {1 — T)K which accounts for this. Our final parametrix K has 
the form 

k ={1- T)K + RtJ 
where the operator Rt verifies the following properties: 

(79) ||V-RT/|||D,|''i- = LPiL9inX'' + II-Rt/IIxs+1 < 11/11 iz)^i-P2-lP2L52 

and 

(80) \\{PRT-T)f\\ ^ \\f\\\D4-P-^Lp' Li' ■ 
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For Tf we have the representation 

Tf = J2Mx)fk{t), hit) = Skf{t,0). 

Then we define 

fcez j>k 

Here S** is a Littlewood-Paley decomposition in the time-frequency variable. That is, 

oo 

1 = E ^KA) 

j=-oo 

with 

supp 4 C {2^-^ < |r| < 2-'+^}. 

The notions of Si,^^ S'<j, etc. are then analogous to those defined in Section O denotes 
the operator with Fourier multiplier r~^, where r is the frequency variable dual to t. 
For fk we estimate 

Since 

n n 1 

-P2 + — = 77 + — , 

we obtain 

(81) IIML.^^2(-^+t+^)l5,/|||^^l_,^_.^,,^,,. 

We now proceed to estimate Rrf ■ By Bernstein's inequality in time we have 

||V(,#.,+i(i) -^j(i))Dr'S>,A(i)ll|D,|.,-.i.iM < 2-iT'2-'2l-«M2'4"*"||A||^,. 

= 2-'-t+i"||A||^., 

S2'-'*»*A'"-"||S./||,„_,_„_.,,.,.,. 

Since 

n — 1 n 1 

s < < - H , 

2 - 2 p2 

it follows that we have off-diagonal decay, while the diagonal summation is controlled by the 
Littlewood-Paley theory. This works if gi 7^ cxd. In the special case gi = 00 we also need 
to observe that the bump functions — concentrate in different spatial regions; 

therefore cannot produce pointwise accumulation. 
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We continue with the norm: 



jZ)^|-P2-sLP2L«2 • 

and the summation works out as before. 

The X*+^ norm is next. Taking advantage of the fact that — 4>j){0) = we compute 

where the restriction s < insures that the norm on the left is finite. We still have 
off-diagonal decay, and for the diagonal summation we can use spatial orthogonality. This 
concludes the proof of ( 1791) . 

For the error estimate (IHOj) we split 

P = -Dl + A. 

The expression PiRtJ is bounded in the same manner as above. On the other hand we have 

k j>k 

and for the summand on the right we can use again Bernstein's inequality with respect to 
t. □ 



Now we prove (1231) . If 

Pau = f + g, fe ID.r^^-^P^L^^ g e r^ 

then we write 

u = Kf + V. 

We use (ITTi) to bound VKf in \Dx\^^^''^LP'^ L'^^ fl X^. It remains to bound v, which solves 

PaV={l-PaK)f + g. 

In the case of Theorem [6] we use successively f l73l) . Theorem HI flTTj) . and fl72l) . We obtain 

< ||Vt;(0)||^. + ||Pat;||y. 

< ||V«(0)||^. + WVKfW^^^. + 11(1 - P,ir)/||y. + \\g\\ys 



< l|VMO)||^. 



Dx\-P2-'LP2L''2 
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This establishes fl23l) with P replaced by Pa- Using f|T9l) . fl23|) then follows. 

In the case of Theorem [7] the argument is similar, but instead of using Theorem H] we 
assume that the localized energy estimates hold. □ 



8. PSEUDODIFFERENTIAL OPERATORS AND PHASE SPACE TRANSFORMS 

Here we tersely introduce the microlocal setup which will be required in the sequel. A 
more detailed exposition can be found in [36], [37j, and the references therein. 

Precisely, our initial goal is to provide a phase-space description of the flow for a pseudo- 
differential evolution of the form 

(82) {Dt + a'^{t,x,D))u = 0, u{0) = uq 

with a real symbol a. We begin by introducing a simpler set-up, which suffices in order to 
obtain a short time description of the flow. In terms of symbol classes, we begin with the 
standard class Sqq of symbols a satisfying 

|a:afa(x,OI <c./3, \a\ + \P\>0. 
We also need the following generalizations S''^'^^ = Sq^''-' of the above class, defined by 

\d:d^^a{x,0\<Ca(3, \a\ + \(3\>k. 
For a phase space transform we use the Bargman transform T defined by 

Tu{x,0 = Cnj e-^e'^^^-y^uiy) dy. 

This is an isometry from L^(]R") to ^(M^") and thus satisfies T*T = I. However, T is not 
an isomorphism; instead, its range consists of functions which satisfy the Cauchy-Riemann 
type equation 

(83) id^T = {d, - ii)T. 

To each pseudodifferential operator a^{x,D) we associate its phase space kernel, i.e. the 
kernel of the conjugated operator Ta^{x, D)T* . A simple example of the correspondence 
between the symbol class and the phase space kernel is the relation (see [36], Theorem 1]) 

a G ^ \K{{x,0, {y,v))\ < cn{i + |(x,0 - {y,vW V iV G N. 

This leads to an easy proof of the Calderon-Vaillancourt theorem, which asserts that the 

operator a"" is L2 bounded if a G S^^\ 

We now turn our attention to the equation (l82l) where we assume a is a real symbol, is in 

S*^^-* uniformly in t G [0, 1], and is continuous in t G [0, 1]. This suffices in order to guarantee 

that (!82l) is well-posed in L^. We let S{t, s) denote the evolution operators corresponding to 
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fl82|) : these are all isometries. We denote by K{t, s) the phase space kernels of S(t, s), i.e. 
the kernels of TS{t, s)T*. It is natural to try to characterize the kernels K{t, s) in terms of 
the Hamilton flow associated to fl82|) : 

^g4) \x = a^{t,x,0 

\i = -a^(t,x,0- 

The corresponding phase space evolution is denoted by x{t, s). These are canonical transfor- 
mations in M^". Furthermore, the condition a E S*^^-* guarantees that s) are bilipschitz 
uniformly with respect to (t,s) G [0,1]. As it turns out, the phase space kernel K{t,s) of 
S{t, s) can indeed be easily characterized as follows: 

Proposition 23. [36l Corollary 7.4] Assume that a is a real symbol in S^'^^ uniformly in 
t G [0, 1]. Then the phase space kernels K{t, s) of S(t, s) satisfy 

\Kit,x,^,s,y,v)\ < C7v(l + \{x,0 - xit,s)iy,v)\r''- 

We also have a corresponding Egorov theorem. For a pdo q^{0), we define its conjugate 
with respect to the flow by 

g'"(t) = S(t,0)g'"(0)5(0,t). 
Then the counterpart of Egorov' s theorem in this setting is 

Proposition 24. [36, Proposition 7.6, Proposition 7.7] Assume that a is a real symbol in 
S^'^^ uniformly in t E [0, 1]. 

(a.) If q{0) E S^^\ then q{t) E 5*^°^ uniformly in t. 

(b.) If q{0) E S^^\ then q{t) E S^^"* uniformly in t, and 

g(t,x,0-g(0)ox(0,t)G5(°). 

The counterpart of this result for g(0) G S^"^^ is not valid in general. However, we can 
prove it in a special case, which will be useful later. 

Proposition 25. Let A > 1. Assume that a{t,x,^) = \\^\, and let q{0) E S^"^^ be an operator 
which is localized at frequency A. Then, q{t) E S^"^^ uniformly in t E [0, 1] and 

g(t,a;,0-g(0)ox(0,t)G5(°). 

Here, analogous to the definition in Section 2, we say that an operator K is localized at 
frequency A if both J{f and K*f are supported in {2-^°A < |^| < 2^°A} for all / G 5(R"). 

We remark that, in the context of the Schrodinger equation, a similar result was proved 
in [37] for a(0 = e'- 
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Proof of Proposition \2B We explicitly compute 

Then we want to show that 

r"'(t, X, D) = e-'*^l^lg"'(0, x, D)e''^\^\ - q'"{x - \tD\D\-\ D) e OPS^^^ 
uniformly in t G [0, 1]. Compute 



where 

rns,x,D) = -tX[\D\,q'^{x- sXD\D\-\D)] - ^q"" {x - s\D\D\-\ D) . 



Using the Weyl calculus, as |,^| ~ A we get 

By Proposition [24l conjugation by e^*'^*'^' leaves the S^^^ class unchanged from which the 
conclusion follows. □ 

From the perspective of the present work, the main disadvantage of Proposition [23] is that 
it can only be used on a fixed time-scale. Of course, appropriate versions can be obtained 
for other time scales simply by rescaling. For instance, in order to obtain results which are 
valid up to time s we need to replace the Bargman transform T with its rescaled versions 

TiM(t,x,0 = CnS--i j e-^^e'^^''-y^uit,y) dy. 

This is often called the FBI transform. It is still an isometry, and its range consists of 
functions satisfying the rescaled Cauchy-Riemann type equation 



^5) -aTi = (d^ - inxi . 

s f s 



Correspondingly, the symbol classes 5**^^^ are replaced by Ss defined by 

a E := {|9^9fa(s,x,0l < c^/js^, |a| + \p\ > k}. 

These are rescaled versions of the S^'^'^ spaces, and thus, results on S^''^ can easily be trans- 
ferred to these classes. In this context, the decay of phase space kernels would be measured 
with a rescaled distance function 
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Still, rescaling does not bring us closer to our goal, which is to work on an infinite time 
scale. This difficulty was resolved in [37] by using a time dependent scale to study the 
evolution (|82l). 

9. A LONG TIME PHASE SPACE PARAMETRIX 

In this section, following ^7\, we consider global in time evolutions of the form 

(86) {Dt + a'"{t,x,D) -tb'^{t,x,D) + c'"{t,x,D))u = 0, t>0 

with time dependent scales for the symbols a, b, c. Precisely, we introduce the classes l^S^'^^ 
of symbols in M x T*M" whose seminorms are given by 

Y^2^{'-^^)\\d:d^^a{t,x,0\\L^i{t^2.}), \a\ + m > k. 
j 

When k = 2 we also need to better track the second derivatives of the symbols using 
the function e{t) introduced in Section [21 We denote by /^Sp^ the subset of I'^S'^'^^ whose 
seminorms are 0(e) when |a| + \f3\ = 2. This additional condition can be rewritten as 

(87) \0-d^^a{t,x,0\<'-^t-^, \a\ + \(3\ = 2. 

We consider the equation (l86l) with a real symbol a G l^Se'^\ which drives the evolution, 
b e l^S^^^ with b>0, which is a damping term, and a possibly complex symbol c G l^S^^\ 
which can be regarded bounded error. 

We remark that the symbols a, b, c above are not related to the coefficients a, b, c of P, 
though they play somewhat similar roles. This slight abuse of notation is harmless since 
at this stage our arguments no longer involve the coefficients a, b, c of P; instead all the 
analysis in the parametrix construction is done on the half- wave evolutions at frequency 1, 
using the symbols a^. 

We let S(t,s) now denote the evolution operator corresponding to (1861) . The following 
result on the evolution was shown in [37] and follows from fairly standard energy estimate 
techniques. 

Proposition 26. [37l Proposition 28] Assume that a G l^S^"^^ and b G l^S^^'^ are real symbols 
with b > 0, while c G l^S^^\ Then the equation (!86|) is forward well-posed in L^(]R"), and 
the corresponding evolution operators satisfy 

\\S{t,s)\\L2^L^<l, 0<S<t. 
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The evolution fl86l) is considered in [37] using a time-dependent phase space transform. 
Precisely, at time t one uses the FBI transform Ti . Thus, the phase space kernels K{t, s) of 
S{t, s) are defined to be the kernels of the conjugated operators 

S{t,s) =TiS{t,s)Tt. 

i s 

A main result in [37] is to establish precise bounds on the phase space kernels K{t, s). These 
bounds are described in terms of the Hamilton flow dictated by a and the decay dictated by 
the damping b. 

The Hamilton flow of Dt + is given by (|84|) and as above, we use s) to denote the 
evolution operators. We shall use 

to denote the trajectories of the flow. Using the linearized equations, one can compute the 
Lipschitz regularity of this flow. See [STJ Proposition 29]. It turns out, however, to be more 
convenient to parametrize s) using the variables [xg, C,t)- In this context, one obtains the 
following regularity. 

Proposition 27. [371 Equation (73)] If a G l-^S), with e small and s < t then 
(gg^ dixt,Q _ fln + eO{l) eO{t) 



In order to describe the decay caused by the damping, we define 

'4'{t,xt,^t) = / b{s,Xs,Q ds. 



We expect b to cause the energy to decay like e '^(*'^'''^*) along the flow. Using the linearized 
flow, it can be shown that 

Proposition 28. [37] Proposition 30] //a G l^S^^^ with e small, b G l^S^^^ and t > s then 

(89) aw...i.)-,p(..,Q) ^ 1 , 



In terms of the above quantities, we can now state the pointwise bound on the kernel of 
the phase space operator S(t, s). This is one of the principal results of [37] . 
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Theorem 29. [371 Theorem 31] Let a G l^S^^\ b G l^S^^^ be real symbols with b > and 

c G l^S^^\ Then for s < t the kernels K{t, s) of the operators S{t, s) satisfy the bound 

/ I _ |2 ^ 

(90) \K{t, X, 6, s, Xs, 1 <t-^-sUl + a - ij{xt, 6))' + + s\^-^s 



To prove this result one considers the phase space evolution of v{t) = Tiu{t) where u solves 
( !86|) . As it turns out, modulo negligible errors this evolution is governed by a degenerate 
parabolic equation with the following components: 

(a) A transport term along the Hamilton flow of a 

(b) A damping term produced by b 

(c) A degenerate parabolic term which is due to the change of scale in the FBI transform. 
Pointwise bounds for the kernel of this degenerate diffusion are obtained in [37] using the 

maximum principle. 

10. A PERTURBATION OF THE HALF WAVE EQUATION 

The results in Theorem apply for symbols a which satisfy the smallness condition 
a G l'^Se'^\ Instead, the symbols small perturbation of ±|^|, precisely 

a± G ±1^1 ^ 1, \x\ ^t. 

To remedy this, in this section we consider the evolution equation 

(A + \D.,\ + a^{t, X, D) - ib'^it, X, D))u = 

where Oq G l^sf\ bo G l^S^^^ are real symbols with bo > 0. Since we are interested in 
this evolution only at frequency 1, we will also make the symplifying assumption that Oq 
is localized at frequency 1 and that b^ — bo{t) is also localized at frequency 1. Here bo is 
simply a function of t. These assumptions guarantee that if the initial data u{to) is localized 
at frequencies {|^| G [2~^'^, 2^'^]} then the solution u inherits this localization. The above 
evolution will serve as the model for our outgoing parametrix. 

We denote by So(t, s) the evolution generated by the above equation. Due to the above 
frequency localization of and b^ — b we have 

So(t, s)S'_io< . <io = 'S'_io< . <ioSo(t, s) = S'_io< . <ioSo(t, s)S'_io< ■ <io- 

We denote by 5*0 (t, s) its (frequency localized) phase space image 

So(t,s) = TiSoit, s)S^io<:.<:ioTt- 

t s 
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We want to obtain bounds on the kernel of So{t, s) which are similar to the ones in Theo- 
rem [29l As a preliminary step we need to study the regularity of the associated Hamilton 
flow which we denote by Xo(^) s)- This can be done directly, but for our purposes it is more 
convenient to reduce it to the case considered in the previous section. 
At each time t we consider the symplectic map ^ defined by 

which corresponds to the Hamilton flow for the Df + \D\ evolution. This extends to a 
space-time symplectic map 

fi{t,T,X,0 = {t,T-\C\,X + tm~\0- 

If po is the symbol 

Po{t, T,X,^) = T+ 1^1 + ao(t, X, 0, 

then its image through fi is 

Hence the conjugate of the Hamilton flow Xoit, s) for r + |,^| + ao with respect to fit is the 
Hamilton flow s) for r + a{t, x, 

We note that a G /^S'P'' iff G l^Se'^\ Hence from (IHHj) we obtain its counterpart for the 
Xo flow, 

Proposition 30. If aQ G l^S^^^ with e sufficiently small and t > s then the Hamilton flow 
Xo(t,s) has the Lipschitz regularity 

ran _ f 4 + e0(l) 2(t- + eO(t) \ 



We proceed in a similar manner with bo and set 

b{t,x,0 = bo{t,x + t^\^\-\0- 
Then the integral ipo of along the Xo flow is the fi conjugate of the integral if) oi h along 
the X flow. Hence we also trivially obtain the analog of Proposition [28l namely 

Proposition 31. If G l^S^f^ with e sufficiently small and bo G l^S^^^ then for t > s we 
have 

(92) diM^t,^j -M^s,Q) ^ 

(y{Xs, qt) 
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Now we can state our main result: 

Theorem 32. Let qq ^ l^S^^\ bo G l^S^'^^ be real symbols with bo > with e sufficiently 
small, so that and bo — bo{t) are localized at frequency 1. Then for s < t the kernel Kq of 
the operator So{t, s) satisfies the bound 

(93) 

/ I 12 \ 

\Ko{t, X, 6, s, x„ 01 ^ ^"^s? f 1 + (M^s, Q - M^t, ^t)f + + - 

x(l + tci(|e.|,[2-^2i°]))-^. 



Proof. We use Theorem [22] via a conjugation with respect to the flat half- wave flow, which 
corresponds to the canonical transformations fif Denote 

S{t,s) = e''^^^So{t, s)e-'^l^l. 

Then we compute 

j^S{t,s) = -ie''\^\{-a^{t,x, D) + ib^{t,x, D))e-''\^\S{t, s). 

Hence S(t, s) is the evolution associated to the pseudodifferential operator 

e''^^^{-a^{t, X, D) + ib'^it, X, D))e-'*^^^. 

Using rescaled versions of Propositions [2^11251 this operator can be expressed in the form 

a"'(t, X, D) - ib'"{t, X, D) + c"'(t, X, D) 

where the remainder term satisfles c G l^S^^\ Hence the phase space kernel of ^(t, s) satisfles 
the bounds given by Theorem [291 

Returning to the original equation, for the phase space evolution S'o(t, s) we can write 
5o(t, s) = Ti e-^*l^l 5_io< . <io5(t, s)5_io< . <ioe*'l^'Tr 

t t t s s s 

= (Tie-**I^IS_io<.<ioTr)5(t,s)(Ti5_io<.<ioe^'l''lTr). 

t t s s 

By a rescaled version of Proposition [231 the kernel of the flrst factor Tie"**'^'S'_io< . <iori is 

11 * t 

rapidly decreasing on the t^ x t~2 scale away from the graph of fit as well as away from the 

support of the symbol 5'_io<.<io, while the kernel of the last factor Tie*'^l^l5'_io< . <io7'r is 

11 s s 

rapidly decreasing on the 52 x s~2 scale away from the graph of fi~^ as well as away from 
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the support of the symbol 5'_io<.<io- Hence the composition simply replaces the Hamilton 
flow associated to a by the Hamilton flow associated to Oq and the function with ipQ in the 
kernel bounds. Thus fl90|) implies (l93l) . and the proof is concluded. □ 



11. The parametrix construction 

We end with a proof of Proposition[T51 That is, we construct a parametrix Kq for Dt+A^^^y 
and easy modifications yield also a parametrix for Dt + A^-j . In the sequel, we shall drop 
the + signs and denote these by Kq and Dt + A{q) respectively. The ± signs will be reserved 
to distinguish waves which are outgoing forward, respectively backward, in time. 

We partition the annulus |^| ~ 1 in phase space 

± i>o 

with 

supp pf C {2^-1 < \x\ < 2^+\ ±x^ > -2-^\x\}, j > 1, 

supp Pq C {|x| < 2, ±x^ > -2^^\x\}. 

At the expense of Schwartz tails which play no role in the sequel, we may replace p^ by 
S^^io{Dx)pf. As such, we shall do so without changing the notation. This allows us to 
assume that the operators are frequency localized to frequency 1. 

In the proposition which follows, we construct evolution operators S^{t,s) as the evolu- 
tions associated to a certain damped half-wave equation. We then form Kq by setting 

^u.^fE';L,S-it,s)iPrnx,D), t<s 

\ET=isnt^^)iprris,D), t>s. 

The properties of Kq listed in Proposition [15] follow easily, after summing, from the given 
properties of S^. 

Proposition 33. Assume that e is sufficiently small. Then for each s G M, there is an 
outgoing parametrix for Dt + A{q) in {t > s} which is localized at frequency 1 and 
satisfies the following: 

(i) bound: 

\\Sl{t,s)U2^L^<l 

(ii) Error estimate: 



(94) ^ ^ 



-N 



x'^Dt{Dt + A(o))5+(t, < (2^- + \t- sl)-^ 
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(iii) Initial data: 

S+{s + 0,s) = I 

(iv) Outgoing parametrix: 

(95) ||l{|,|<2-io(|i_,|+2.)}5';(t, s)P;\\l2^l2 < {\t -s\+ 2^)-'' 

(v) Finite speed: 

(96) ||a:n{|,|>2io(|t-.|+2.)}5;(t, s)P;\\l^^l2 < {\t - s\ + 2^)-'' 

(vi) Frequency localization: 

(97) 11(1 - PM,4])^;(i, s)P;h2^L'^ <i\t-s\+ 2^^" 

(vii) Pointwise decay: 

(98) ||5'+(t, s)P;\\li^l^ < (1 + |t - sl)-"^. 
With obvious modifications, the same hold for Sj'. 

By translation invariance, without any loss of generality we may assume that s = 2^ . We 
first reduce the problem to the study of an evolution of a perturbed half-wave equation as 
in Section [101 Heuristically we observe that in the support of the symbol of P^ we have 

G [2-2,2^], \x\^s, x-^>-^\xm. 
An easy computation shows that along the forward Hamilton flow starting here we have 

1^1 G [2-\2% \x\^t. 

But in this region we have 

a(o)(t,x,0-|e|e/iSf) 

which follows from the analog of (132!) which holds for a(o). Thanks to (1951) . (196!) and (197|) . 
we can freely modify the symbol of a(o) in the regions <^ t} and {|,^| ^ [2~^,2^]} at 
the expense of producing a negligible error in ( l94l) . It thus suffices to study the evolution 
governed by a symbol 

|e|+ao(t,x,0, aoG/^^f 
so that flo vanishes if {|^| ^ [2~^,2^]} and is localized at frequency 1. 

In essence, = a(o) — |^|, and thus by (1521) . we may assume the better regularity 

|a:afao(t,x,OI <e(t)rl°l, |a| <2 



(99) 



\d:d^^ao{t,x,0\<e{t)t-'-^, \a\>2. 
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This additional decay shall be used on time scales which are too small to allow x 
packets at time s to separate in time t. 

Unfortunately, simply defining the parametrix by the evolution associated to the op- 
erator \D \ + 0(0) does not seem to work. Precisely, the bounds (!95|) . (!96l) and (!97|) appear 
to fail. This is because at each time t, there is leakage caused by the uncertainty principle 
to the regions appearing in ([HSD, and fp7|) . which are outside the propagation region 
indicated by the Hamilton flow. While this leakage does have rapid spatial decay, its time 
evolution yields output which does not have the rapid decay in time as needed in ( 195|) . (I96l) 
and (EI]). 

Thus, in order to be able to prove the rapid t-decay in, e.g., (j95l) (!96l) and (!97l) . we shall 
introduce an artificial damping term bo E l^S^^\ bo{t,x,^) > 0. The role of 60 is precisely to 
put a damping on the time evolution of the above mentioned leakeage. At the same time, 60 
is taken to be in the main propagation region. We would hke to define Sj'{t, s) to be the 
forward evolution operator of the equation 

(A + 1^1 + a^{t, X, D))u = tb^{t, X, D)u. 

However, in order to insure the frequency localization of our parametrix we replace s) 
by the truncated operator 

We shall show that 

(100) ||x-5<_5(A)5+(t, s)PI\\l2^l^ < {\t -s\+ 2^)-'', 

(101) \\x'^d^S^,iD^)Slit, s)P;\\l2^l^ < i\t -s\+ 2^-)-'^, 

and thus, the errors in flMl) which result from this truncation are negligible. We shall further 
prove the following bound on the damping term 

||x"6^(t, X, D)S^{t, s)P^\\l2^v^ < (|t - + 2^r'^, 
102 11 

Wx'^Dtb'^it, X, D)S^{t, s)P;\\l2^l^ <i\t-s\+ 2^)-^, 
which shall yield flMl) . 

With Sj'{t,s) now fixed, property (iii) is trivial, and (i) follows from Proposition [261 We 
proceed to the argument which yields our main pointwise bound (!98l) . Here, we examine 
three cases separately. 

Case 1: |t — s| > s. In this regime, we may neglect the damping. For initial data 
u{s) = 60, we have 

^2 

Tiu{s,Xs,0 = s~?e"2fe""^. 

s 
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Using Theorem [32], we see that 

|TiM(t, X, 6) I < r ^ y"(l + t-^\x - xt{it, X,) + s\i - U^s, it)?r''e-^^ dxs 

<t-?s-t j {l + t-^\x-xt{iuxs)\^ye-^^dxs. 

For the remaining integral, we use that Xs Xt{C,t,Xs) is Lipschitz. See ( lOTll . Integrating in 
Xs then yields 

a;, et)| < r?(i + r V - o)|2)-^, 

and by applying Tt, we have 

Ht,y)\<t~^ j (i + r'\x-xt{^t, 0)1^)-'' dxd^t 
< /(i + ri|^-x,(6,o)p)-^t/6. 



If \t — s\ > s, then the map C,t Xt{C,t, 0) is zero homogeneous, by ( 19T1) has Lipschitz constant 
which is bounded by t, and has maximal rank n — 1. Hence, integration with respect to 
yields 

Ht,y)\<t-'^. 

Case 2: 1 < \t — s\ < s. Here, we reinitialize the time scale to prevent difficulties which 
result from x s"^ packets at time s not separating before time t. In addition to (199!) . we 
similarly require 

|9:9f6o(t,a:,OI <t-^-'"', I«I<1, 



\d:dX{t,x,o\<t-'-'^, i«i>i 



|c';;fc'^'6n(t,x,e)| < t-' 

(103) 

for |t — s| < 2-^. The additional regularity (p9l) and fll03p is sufficient to show that Oq, &o 
remain in the appropriate symbol classes after the time translation which sets the initial 
time to t — s. Theorem [32] thus remains valid, and the bound follows from the computation 
above in the translated coordinates. 

Case 3: < |t — s| < 1. Here, since our initial data is localized at frequency 1, we may 
simply use Sobolev embeddings combined with the bounds from Proposition [221 

||5+(t,s)P/no||Loo < \\S^{t,s)P;uo\\v^ < ||P>o||l2 < holUi. 

The rest of the proof is based on properties of Bq. In particular, we use a construction 
which is quite similar to that of [37] to build a Bq which allows us to prove the remaining 
required estimates: (195]) . (196]) . fllOOl) . fllOip . and fll02p . In particular, we have 

Lemma 34. There exists a symbol b E l^S''^^ which satisfies, in addition to fll03p . 

48 



(bl) t^b is nonincreasing along the Hamiltonian flow for Dt + \Dx \ + , and 

< th{t,Xt,it) < 1 =^ b{2t,X2t,^2t) = 0. 

(b2) At the initial time, we have 

6(2^x,0 = 0, zn {2-3 < |^| < 2^, 2^-^ < |a;| < 2^+\ > -2-^\x\}. 
(b3) At any time t > 2^ , we have 

b{t, X, = t''^, outside {2"^ < |^| < 2^} n {2-^^ < \x\ < 2H}. 

Before proving this lemma, let us explain how such a damping term can be used to complete 
the proof of Proposition [T51 Indeed, we have the following lemma which is essentially from 

Lemma 35. Assume that the symbol bo E l^S^^^ satisfies the properties (bl), (b2), and (b3) 
above. Then, the bounds ([95]), ([96]), ([97]), (fTOOD . (fTOT]) . and ffT02D hold. 

Indeed, once Theorem [32] has been established, the necessary modifications to the argu- 
ments from [3^ are quite simple. It only remains to establish the second estimate from (11021) 
which did not appear in [37]. Here, however, we note that, modulo negligible errors due to 
the frequency truncation of S*"*", 

A&o a;, D)S+{t, s) = -tidAnt, X, D)S+{t, s) - b^{t, x, D)\D,\S+{t, s) 

- b^{t, X, D)a^{t, X, D)S+{t, s) + ib^{t, x, D)b^{t, x, D)S+{t, s). 

Since the symbols of (dtbo)"', b^lD^l, &o '^o ' ^o^o ^t^'^ have supports which 

are contained in the support of bo, we may similarly apply Proposition 17 of [37] to obtain 
the estimate. 

It now only remains to complete the construction of said damping terms b. 
Proof of Lemma^S^ We define the increasing bounded function e{s) by 

e(.) = r '-^da. 
Jo o- 

Letting be a smooth, nondecreasing cutoff function which equals in (— oo, 0) and 1 in 

(1, oo), we set 

b{t, X, = r i (1 - 0(6l)0(&2)0(&3)0(&4)0(&5)) 

with 
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• Cutoff frequencies which are too large 

b(tn- ^'^' + <t)-\^\ 

• Cutoff frequencies which are too small 

\C\-2-V^ + ce{t) 

7(t) 

where c is a fixed small constant, 

• Select outgoing waves 

• Cutoff values of \x\ which are too large 

2H - \x\ 



2-^1 




2-i2| 


\x\ 





t 

• Cutoff values of \x\ which are too small 

mi - 2-H\^\ + x^ 

We note that 



{2-=^ < 1^1 < 2^} n {2-H < \x\ < 2H} n {x^ > -2-^]} = A c {6 = 0} 

if e is sufficiently small, while 

{th <i}cEt = {2^4 < lei < 2^} n {2-H < \x\ < 2H} n {x^ > -2~'/^\x\\^\}. 

So, the conditions (b2) and (b3) arc easily satisfied. 

To prove (bl), it suffices to study the behavior of b along the Hamilton flow within Et and 
show that for each bj, we have 

(104) j^b,{t, xt, 6) > p in Etr){0< bj < 1}. 

Here t — > {xt,it) now denotes a trajectory of the flow for Dt + \Dx\ + a^. For {xt,^t) G Et, 
we have 

We simply calculate 

for e sufficiently small. The computation for 62 is identical. For 63, we have 

d, , .^ \^t\^\xt\'' - (xS? 0(e(t)) 2 ^ , 
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For 64, we compute 



d 



)5 



—bJt, xt) = ' ' ' ' , + ' ' > — , in < 64 < 1- 
Finally, for 65 we also compute 

> ^ > J, in n {0 < 65 < 1}. 
It remains to verify fllOSp . and hence that b E l^S^^\ but this is straightforward. □ 
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